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ABSTRACT 

The  multi-institutional  interdisciplinary  research  program  of  the  Steel  Research  Group  (SRG)  has  centered  on 
the  development  of  a  systems  approach  for  the  design  of  advanced  high-performance  steels  using  computer-aid^ 
modelling.  Initial  long  range  studies  involved  the  development  of  the  theoretical  foundations  underl5dng  strength, 
toughness,  and  hydrogen  resistance.  Mechanistic  based  computational  models  include  the  prediction  of  martensite 
start  temperatures  in  high  alloy  steels,  multicomponent  coarsening  kinetics  in  non-ideal  non-dilute  systems,  and 
the  martensitic  transformation  volume  change  in  multicomponent  systems.  Modelling  has  successfully  related  the 
peak  hardness  in  secondary  hardening  steels  to  the  coherent  ^ving  force  of  the  M2C  carbide  in  para-equilibnum 
with  transient  cementite.  Additionally,  the  precipitation  time  scale  has  been  related  to  the  multicomponent 
coarsening  coefficient  of  the  M2C  carbide  with  respect  to  the  half-completion  state.  Within  a  systerns  design 
framework,  these  models  have  been  successfully  applied  to  the  design  of  high  performance  alloys  including 
stainless  bearing  steels,  advanced  armor  steels,  and  high-hardness  gear  steels. 
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INTRODUCTION 


As  an  innovator  in  the  use  of  computational  metallurgy,  Prof.  Jack  Kirkaldy,  was  an  early  pioneeyn  the  use  of 
computer  power  applied  to  both  the  fundamental  study  and  commercial  development  of  steel.  His  work  has 
inspked  a  generation  of  researchers  to  dedicate  their  careers  to  this  field.  Recendy,  Ae  ^mputational  modelling  of 
phase  transformations  occurring  during  thermal/mechanical  processing  of  steels  has  become  more  complex  ye 
less  expensive  as  computational  capabilities  have  increased  and  the  cost  of  computing  power  has  deceased  At 
the  same  time,  experimental  investigative  techniques  have  also  advanced,  utilizing  the  latest  developtnents  in 
characterization  methods,  albeit  at  ever  increasing  cost.  This  new  relative  cost  structure  ^mpumtio 

modelling  and  empiricism  is  driving  a  revolution  in  the  way  m  which  steels  ^e  developed.  Historically,  me 
development  of  steds  foUowed  a  largely  empirical  path  in  which  quaUtative  models  drove  empincd  stupes  lea^g 
to  the  development  of  a  new  alloy.  The  new  paradigm  of  alloy  design  reverses  these  trachtional  roles. 
Computational  modelling  can  now  perform  the  detailed  simulations  of  the  inicrostructural  res^nse  to 
thermal/mechanical  processing  required  to  determine  a  new  alloy  composition  while  empmcism  vali^tes  the 
computer  models  and  determines  the  optimum  processing  of  a  new  alloy.  Thus  empirical  studies  now  ^ve  Ae 
development  of  computational  models  leading  to  the  design  of  a  new  steel.  Ferrous  naetallurgy,  containing  a 
wealth  of  prior  mechanistic  and  empirical  knowledge,  possesses  a  distinct  advantage  for  the  development  ot 

computational  modelling  and  is  at  the  forefront  of  this  revolution.  ,  .  j  t  *u-  u 

Development  of  the  computational  models  must  reflect  the  manner  in  which  they  are  to  be  used.  In  tois  work, 
the  modelling  is  part  of  a  general  systems  approach  to  the  design  of  materials.  The  systems  design  of  materials 
builds  on  the  universal  systems  engineering  approach  described  by  G.  M.  Jenkins[l]  ^d  the  concept  of  matenals 
as  systems  by  Cyril  Smith[2]  describing  material  structure  as  a  hierarchy  of  interacting  microstructu^ 
subsystems.  This  approach  is  further  expanded  to  describe  the  hierarchial  nature  of  processing  and  materim 
properties[3]  using  a  linear  processing-structure-properties-performance  model  of  maten^s  science  depict^  by 
the  three-link  chain  model  of  Figure  1.  An  example  of  this  Merarchial  descnption  is  given  in  the  Aow-bltxk 
diagram  of  secondary  hardening  ultrahigh-strength  steel  of  Figure  2.  The  flow  block  dia^ain  illustrates  the 
influence  processing  steps  have  on  the  microstructural  components  of  the  matnx,  stiengthenmg  dispereion, 
austenite  dispersion,  grain  refining  particles  and  grain  boundary  chemistry  and  the  resulting  desired  properties  ot 
strength,  toughness  and  hydrogen  resistance.  Describing  each  of  these  components  as  a  set  of  mteractmg 
subsystems  arranged  in  a  hierarchial  manner,  computational  modelling  serves  the  role  of  providing  quantitative 
descriptions  of  the  interactions  between  subsystems.  In  order  to  be  successful,  the  modelling  must  be  mech^istic 
and  quantitative  in  nature.  Mechanistic  computational  models  can  be  used  as  predictive  design  twls,  capable  ot 
extrapolation,  while  empirically  based  models  are  interpolative  and  limited  to  problems  in  which  the  solution  is 
known  to  lie  within  the  bounds  of  the  model. 


Figure  1 


Thiee-link  chain  model  of  materials  science  illustrating  the  linear  processing- 
stnicture-properties-performance  description. 


PROCESSING 


STRUCTURE 


PROPERTIES 


The  Steel  Research  Group  (SRG),  a  multi-institutional  research  program  compnsed  of  members  from 
industry,  academia,  and  government,  was  initiated  in  1985  to  combine  physical  metallurgy,  applied  mechamcs, 
physics,  and  chemistry  in  order  to  establish  the  basic  principles  underlying  the  processing-stnicture-proi^rties 
relations  in  steels  and  apply  this  understanding  to  the  design  of  superior  alloys.  The  goal  of  this  effort  consists  of 
the  development  of  broad  based  quantitative  relationships  to  describe  the  processing-structure  and 
structure-property  interplay  and  a  systems  approach  to  using  these  tools  that  lead  to  the  CTeaUon  of  unprov^ 
alloys.  The  development  of  quantitative  models  for  the  prediction  of  alloy  behavior  is  well  under  way  and  the 
recent  focus  is  centered  on  the  framework  in  which  to  use  these  tools  to  improve  alloys.  As  this  framework  is 
developed,  new  insights  into  the  basis  of  the  modeling  tools  are  found  and,  subsequently,  me  models  ^ 
improv^.  In  this  manner,  a  coherent  approach  to  the  understanding  of  given  matenal  systems  is  achieved  that 
allows  direct  application  to  technologically  important  problems.  Recently,  the  Ferriuin  Division  of  Northwestern 
University’s  BIRL  Industrial  Research  Laboratory  was  established  to  pursue  the  apphcation  of  the^  mo^ls  ^d 
techniques  to  the  development  of  proprietary  alloys  and  processes  for  industry.  This  paper  will  review  key 
computational  modelling  efforts  of  the  SRG/Ferrium  Division  and  provide  an  example  of  how  they  are  used  in 
alloy  design  using  the  systems  approach. 

COMPUTATIONAL  MODELS 


A  vital  component  of  the  modelling  we  have  undertaken  is  the  ThermoCalc  thermochemical  databank  system(4) 
developed  at  the  Royal  Institute  in  Stockholm.  The  Them^alc  program  uses  therm^ynamic  assessments  tom 
binary,  ternary,  and  quaternary  systems  to  extrapolate  to  higher  order  systems.  Equilibria,  constrained  eqmlibna, 
and  driving  forces  can  be  calculated  as  functions  of  composition,  chemical  potential,  as  well  as  user^  denned 
functions.  Central  to  the  versatility  of  the  program  is  a  separation  of  the  calculation  package  from  the  datable, 
simplifying  the  ability  to  use  different  thermochemical  models  to  describe  the  thermodynamics  of  the  system.  TTie 
ThermoCalc  software  offers  unprecedented  access  to  equilibrium  thermodyn^ic  quantities  for  complex, 
high-order  systems.  To  apply  this  information  to  the  modelfing  of  highly  nonequilibrium  processes  of  interest  in 
real  alloys,  we  describe  the  dynamic  nature  of  phase  transformations  in  terms  of  thermodynamic  scaling  factors 
which  can  then  be  evaluated  by  the  software.  The  first  computational  model  that  we  review  is  based  on  this 
approach  and  involves  the  precipitation  of  M2C  carbides  leading  to  the  secondary  h^dening  response  in 
ultrahigh-strength  steels.  The  second  modelling  effort  we  shall  review  involves  the  application  of  first  principles 
quantum  mechanical  calculations  to  the  design  of  steels  with  improved  resistance  to  hydrogen  embrittlement  and 
intergranular  fracture. 


Spcondarv  hardening 

Ultrahigh-strength  (UHS)  secondary  hardening  steels  are  strengthened  by  the  precipitation  of  coherent  MjC 
carbides  during  tempering.  In  high  Co  steels  in  which  dislocation  recovery  is  retarded,  the  M2C  carbides 
precipitate  coherently  on  dislocations  and  provide  the  characteristic  secondary  hardening  peak  dunng  tempenng. 
The  SRG  conducted  a  deliberate  and  comprehensive  investigation  of  this  phenomenon  in  commercial  and  model 
alloys  in  order  to  gain  a  fundamental  understanding  of  the  prixess  and  develop  a  quantitative  design  mc^el.  A 
wide  range  of  techniques  were  utilized  to  gather  experimental  information  across  the  complete  range  of  size  and 
time  scales  of  interest.  Atom-probe  field-ion  microscopy  (APFIM)[5,6],  transmission  electron  microscopy 
(TEM)[7],  small  angle  neutron  scattering  (SANS)[8],  and  X-ray  diffraction  (X^)[7,9]  techniques  provided 
information  on  the  size,  shape,  composition,  and  lattice  parameters  of  the  M2C  precipitates  as  well  as  tlw  resulting 
hardness  values  spanning  tempering  times  of  less  than  an  hour  to  more  than  a  thousand  hours.  This 
identified  that  the  precipitation  was  well  described  by  the  theory  developed  by  Langer  and  Schwartz[10]  for 
precipitation  at  high  supersaturation  in  which  the  growth  regime  is  suppressed  and  precipitation  occurs  by  a 
process  of  nucleation  and  coarsening,  maintaining  a  particle  size  close  to  the  critical  size.  . 

Based  on  these  investigations,  two  important  scaling  factors  are  identified.  The  initial  critical  nucleus 
determines  the  size  scale  of  the  precipitates  throughout  the  precipitation  reaction  and  the  coarsening  rate  constant 
dctcnnincs  the  precipitation  time  scale.  Both  of  these  factors  are  vital  to  the  successful  design  of  a  UHS 
The  peak  hardness  in  a  UHS  steel  commonly  occurs  at  the  particle  size  corresponding  to  the  transition  from 
particle  shearing  to  Orowan  bypass.  It  is  also  advantageous  to  bring  the  M2C  precipitation  to  completion  in  order 
to  dissolve  all  of  the  transient  cementite  which  otherwise  limits  toughness.  Therefore,  the  smaller  the  initial  critics 
particle  size,  the  closer  completion  of  precipitation  occurs  to  peak  hardness  and  more  efficient  strengthening  is 
obtained.  The  time  scale  of  precipitation  is  also  important  due  to  the  kinetic  competition  between  the  secondary 
hardening  reaction  and  the  segregation  of  impurities  to  the  prior  austenite  grain  boundaries  leading  to  intergranular 

embrittlement  ...  .  . 

The  initial  critical  nucleus  size  scales  inversely  with  the  thermodynamic  driving  force  for  precipitation.  In  the 
case  of  the  M2C  carbide  it  is  important  to  include  the  influence  of  prior  cementite  formation  and  coherency  on  Ais 
quantity.  The  coherency  elastic  self  energy  can  be  evaluated  by  the  calculation  of  an  anisotropic  ellipsoidal 
inclusion  using  the  equivalent  eigenstrain  method[l  1]  and  the  impact  of  solute  redistribution  on  the  resulting  stress 
distribution  is  addressed  by  using  open-systems  elastic  constants[12].  By  relating  the  coherency  strain  to 
composition  via  the  compositional  dependence  of  the  particle  and  matrix  lattice  parameters,  the  compositional 
dependence  of  the  elastic  self  energy  is  determined  in  a  form  compatible  with  the  ThermoCalc  software,  pie 
linear  elastic  self  energy  calculation  represents  an  upper  limit  and  a  correction  factor  is  used  to  fit  the  pi^pimtion 
composition  trajectories  of  a  large  set  of  experimental  alloys.  The  impact  of  prior  cementite  precipitation  is 
accounted  for  by  the  calculation  of  the  coherent  driving  force  in  the  presence  of  the  carbon  potential  due  to 
para-equilibrium  cementite.  This  para-equilibrium  carbon  potential  is  defined  by  m  equilibrium  between  the 
matrix  and  cementite  in  which  the  substitutional  species  are  held  constant  and  only  the  interstitial  carbon  is  allowed 
to  partition.  In  this  approximation  the  cementite  acts  as  a  carbon  source  at  constant  chemical  potential. 

Figure  3  represents  the  level  of  agreement  of  the  strengthening  response  of  the  model  alloys  with  the  above 
model.  The  model  alloys  contain  16  wt%  Co,  5  wt%  Ni,  and  0.24  wt%  C  with  varying  amounts  of  the  carbide 
formers  Q,  Mo,  and,  in  a  few  cases,  V.  The  nickel  content  is  chosen  to  eliminate  austenite  precipitation  during 
tempering  which  otherwise  complicates  the  hardening  response.  In  Figure  3,  the  peak  harness  during  tempering 
at  510  C  is  plotted  against  the  driving  force  for  precipitation  of  the  coherent  M2C  carbide  in  the  presence  of  para- 
equilibrium  cementite.  The  open  circles  represent  alloys  containing  V.  The  relationship  demonstrates  the  ability  to 
predict  peak  hardness  values  within  approximately  +/-  25  VHN  in  this  class  of  alloys.  Work  continues  to 
correlate  driving  force  directly  to  critical  particle  size  and  interparticle  spacings. 

The  time  scale  of  precipitation  at  high  supersaturations,  according  to  the  Langer-Schwartz  treatment,  scales 
with  the  coarsening  rate  of  the  particle  distribution.  The  modelling  pursued  in  this  work  expands  upon  the 
Lifshitz-Slyozov  and  Wagner  (LSW)[13,14]  theory,  describing  the  coarsening  of  spherical  partcles  in  a  binary 
system,  with  the  intent  of  removing  Ae  binary  restrictions  of  the  LSW  theory  and  reforaiulating  it  in  a  manner 
compatible  with  the  multicomponent  thermodynamics  of  the  ThermoCalc  software.  Initially,  this  work  involved 
the  extension  of  the  LSW  treatment  to  ternary  systems[15]  and  was  later  expanded  to  the  more  gen^ 
multicomponent  case[16].  The  result  of  this  analysis  characterizes  the  coar^ning  rate  of  a  particle  of  average  size 
as  a  function  of  the  multicomponent  diffusion  coefficients,  the  equilibrium  partitioning  coefficients,  and  the 
second  derivatives  of  the  Gibbs  free-energy  evaluated  at  the  ^uilibrium  state.  The  surface  energy  and  molar 
volume  are  taken  to  be  composition  independent  and  are  considered  constant.  In  this  form,  the  coarsening  rate 
constant  is  the  result  of  an  asymptotic  analysis  and  is  only  representative  at  very  long  tinte  scales  and  very  close  to 
equilibrium.  This  is  certainly  not  the  case  for  the  precipitation  of  the  M2C  carbide  at  high  supersaturation.  The 
matrix  content  of  alloy  is  far  from  equilibrium  during  much  of  the  precipitation  process,  approaching  equilibriimi 
only  near  completion.  This  effect  is  most  severe  for  alloys  containing  stoichiometric  quantities  of  carbide 
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Figure  3  -  Correlation  of  peak  hardness  during  tempering  at  5 10  C  to  the  coherent  ^ving 

force  for  the  precipitation  of  M2C  carbide  in  the  presence  of  paraequilibrium 
cementite. 


formers  as  measured  by  the  relative  difference  in  the  matrix  aEoy  conteiit  during  precipitation  and  at  equilibrium. 
During  precipitation  in  a  stoichiometric  alloy,  the  alloy  matrix  content  is  of  the  same  order  as  the  overall  alloy 
content,  while  at  equilibrium,  the  matrix  alloy  content  is  very  small.  To  define  a  coarsening  rate  constant  more 
representative  of  the  conditions  present  during  the  precipitation  process,  a  coarsenmg  rate  is  evaluated  at  the  pomt 
when  the  volume  fraction  of  the  precipitate  is  one-half  of  the  equilibrium  value.  This  is  achieved  by  calculating  me 
coherent  equilibrium  for  the  M2C  carbide,  and  then,  adding  energy  to  the  M2C  phase  to  account  for  capillarity 
until  the  amount  of  the  phase  is  half  of  the  equilibrium  value.  The  coarsening  rate  is  then  calculated  tom  the 
thermodynamic  properties  of  this  state.  Figure  4  represents  the  correlation  between  the  precipitation 
half-completion  time  and  the  half-completion  coarsening  rate  constant  of  the  model  alloys  for  which  this  dam  is 
available.  The  experimental  data  spans  a  kinetic  range  of  over  two  decades  and  the  correlation  holds  to  within  a 
quarter  decade  over  the  range. 

Interfacial  Cohesion 

Hydrogen  embrittlement  of  ultrahigh-strength  steels  is  asswiated  with  the  prior  segregation  to  the  grain 
boundary  of  impurities  such  as  P  and  S.  A  thermodynamic  treatment  of  this  phenomenon  by  Rice  and 
Wang[17,18]  illustrates  that  the  potency  of  a  segregating  solute  in  reducing  the  work  required  for  brittle  fracture 
along  a  boundary  is  linearly  related  to  the  difference  in  the  segregation  energy  for  the  solute  at  the  boundary  and  at 
the  free  surface.  Specifically,  a  solute  with  a  higher  segregation  energy  at  the  free  siuface  will  be  an  einbrittler 
while  a  solute  with  a  higher  segregation  energy  at  the  grain  boundary  will  enhance  intergranular  cohesion.  A 
survey  of  reported  segregation  energies  and  embrittling  potency(reported  as  the  shift  in  the  ductile-to-brittle 
transition  temperature  per  atomic  percent  solute  on  the  grain  boundary)  in  Fe-base  alloys  demonstrates  these 
general  trends[17];  however,  the  experimental  difficulty  of  surface  thermodynamic  measurements  gives 
ambiguous  values  for  some  solutes. 


A  study  was  undertaken  by  the  Steel  Research  Group  at  Northwestern  University  to  use  state-of-the-art  total- 
energy  electronic-structure  calculations  and  advanced  supercomputer  processing  to  investigate  the 


electronic  basis  of  intergranular  embrittlement  in  light  of  the  Ri(»-Wang  theory[19].  In  the  study,  f^t  pmciples 
calculations  were  used  to  determine  the  total  energy  of  atomic  cells  representing  the  Fe  X3[110](l  ll)  grain 
boundary  and  ^H)  free  surface  with  a  monolayer  of  an  impurity  solute  present.  The  calculations  were 
accomplished  with  the  full-potential  linearized  plane  wave  (FLAPW)  total  energy  technique.  The  atomic  structi^ 
in  each  case  was  relaxed  to  find  the  minimum  energy  state.  The  results  of  these  calculations  include  not  only  me 
segregation  energies  responsible  for  the  embrittling  or  cohesion  enhancing  effects  of  segregating  solutes,  but  me 
underlying  electronic  structure  of  the  solutes  in  the  boundary  and  surface  environments.  A  comparison  of  the 
directional  covalent  electronic  structure  between  B,  a  strong  cohesion  enhancer,  and  P,  ^  a^  strong  embnttler, 
indicates  the  strong  bonding  of  the  B  atom  across  the  boundary  plane  associated  with  hybridization  of  the  B  2p 
electrons  with  the  Fe  d  band.  This  directional  bonding  is  not  seen  in  the  case  of  the  P  atom  which  does  not 

significantly  hybridize  with  Fe.  .  . 

The  results  of  the  first  principles  calculations  are  correlated  to  the  experimental  embrittling  potency  in  Figure  5. 
The  difference  between  grain  boundary  and  free  surface  segregation  energies,  calculated  by  electronic  structure 
calculations,  and  the  experimentally  observed  shift  in  the  ductile-to-brittle-transition-temperature  are  plotted  for  C, 
B,  P  and  S  solutes.  The  C  and  B  are  shown  as  cohesion  enhancers,  P  and  S  as  embrittlers.  The  computed  energy 
differences  are  in  excellent  agreement  with  the  observed  effects  on  interfacial  cohesion.  T^s  study  is  continuing 
to  investigate  the  embrittling  potency  of  additional  solutes  as  well  as  ternary  solute  interactions  that  may  be  able  to 
alter  the  electronic  stracture  of  embrittling  P  and  S  to  reduce  their  detrimental  impact  on  grain  boundary  cohesion. 

MATERIALS  DESIGN 


Background 

In  order  to  illustrate  the  use  of  such  models  in  the  design  of  advanced  rnaterials  we  will  use  the  example  of 
case-hardenable  secondary  hardening  gear  steels.  Design  considerations  for  high  perfonnance  gears  in  aerospace, 
automotive,  and  other  applications  (hive  gear  technology  to  transmit  more  power  in  less  space  and  weight. 
Current  high  performance  gear  steels  are  typically  stage  I  quench  and  tempered  martensites,  case  hardened  through 


carburizing  to  60  Rc  on  the  surface  with  a  core  hardness  of  typically  35  Re-  Gear  teeth  are  cyclically  subject^  to 
bending  and  contact  stresses.  Failure  modes  in  gears  are  generally  grouped  into  three  cat^ones, 
surface  pitting,  and  subsurface  spalling.  High  surface  hardness  is  used  to  limit  the  twth  bre^age  due  to  bening 
as  well  as  pitting  failures.  If  the  fatigue  strength  drops  below  the  applied  stress 


Figure  5  -  Correlation  between  the  calculated  segregation  free  energy  difference  and  the 

experimentally  observed  embrittling  sensitivity(from  Ref  [17]  with  new  data  from  [19]). 


at  any  point  below  the  surface,  subsurface  spalling  can  result.  Typically  a  1  mm  case  depth  is  required  to  provide 

adequate  fatigue  Strength  to  avoid  subsurface  Spalling.  ^  r^u 

Stage  IV  tempered  secondary  hardening  steels  offer  numerous  advantages  over  conventional  gear  steels.  1  ne 
efficient  strengthening  and  superior  toughness  achieved  in  secondary  hardening  steels  allows  greater  core  and  case 
hardness  to  be  attained  reducing  the  size  and  weight  of  gears  needed  to  transmit  a  given  power.  The  grater 
temperature  resistance  of  secondary  hardening  steels  allows  operation  at  higher  temperatures  and  longer  times 
before  performance  degrades.  In  addition,  recent  results  on  commercial  secondary  hardening  alloys  indicate 
casting  of  these  steels  may  be  possible.  With  these  considerations,  the  Steel  Research  Group  at  Northwestern 
initiated  a  program  to  design  an  advanced  secondary  hardening  case-core  gear  steel  targeted  at  a  70  Rc  s^ace 
hardness  with  a  core  hardness  of  50  to  55  Rc*  With  this  increment  in  material  properties,  a  weight  reduction  w 
gear  systems  of  50%  is  possible  with  new  gear  designs  based  on  these  advanced  steels.  The  efforts  described  in 
this  section  represent  contributions  from  many  sources,  including  undergraduate  design  class  and  senior  projects, 
one  of  which  received  1st  prize  in  the  1994  TMS  Materials  Design  Competition(20). 

Analysis 

The  systems  analysis  of  the  case-core  secondary  hardening  steel  system  is  the  first  step  in  the  design  process. 
Figure  6  illustrates  the  system  structure  for  high  power-density  gears  manufactured  by  three  alternative  processing 
routes,  conventional  forged  ingot  processing,  near  net  shape  casting  and  powder  processing.  Case  hardenable 
secondary  hardening  gear  steels  are  a  subsystem  of  this  flow-block  diagram  and  will  ^  the  focus  of  our 
discussion.  The  sequential  processing  steps  dictate  the  evolution  of  the  case  and  core  microstructures  which 
determine  the  combination  of  properties  required  for  the  overall  performance  of  the  system.  Both  the  case  and 
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Figure  6 


Flow-block  diagram  of  high  power-density  gears  manufactured 
from  case  hardenable  second^  hardening  gear  steel  using  three 
alternative  processing  routes. 


Microsegregation 


core  consist  microstructures  of  a  lath  martensite  with  high  Co  for  dislocation  recovery  resistance  essenu 
efficient  secondary  hardening  and  Ni  for  cleavage  resistance  Strengthening  is  provided  by  the  ^heien 
nrecinitation  of  fine  M2C  carbides  on  dislocations.  This  secondary  hardening  reaction  dissolves  the  tr^sient 
cementite  and  it  is  essential  to  bring  this  reaction  to  completion  in  order  to  eliminate  cemennte  for  high  toughness. 
?h7S  refilg  dispersion  has  a  double  impact  on  toughness.  By  limiting  grain  growth  at  high  temperature 
durii^  solution  treatment,  brittle  intergranular  modes  of  fracture  are  avoided.  The  ^n  refinmg  particles  also  play 
an  important  role  in  the  ductile  microvoid  nucleation  and  coalescence  fracture  ^havior.  Thus  it  is 
prese^e  adequate  volume  fraction  and  size  to  pin  the  grain  boundanes  while  choosing  the  phases  w  th  greatest 
FnSacial  Shesion.  Also  important  is  the  control  of  the  ^ain  bound^  chemistry  to  intergranular 
embrittlement  by  hydrogen  in  association  with  pnor  segregation  of  enjnttling  impunties.  Dunng  tempenng 
impurities  segregate  to  L  grain  boundaries  and  in  the  case  of  P  and  S  reduce  the  mterfacial  cohesion  of  the 
boundary  promoting  intergranular  embrittlement.  A  number  of  methods  can  be  used  to  avoid  this  problem. 
Gettering  compounds  can  be  utiUzed  to  tie  up  the  impurities  in  stable  compounds  r^ucmg  the  segregation  to  the 
erain  boundary.  In  order  to  produce  the  most  stable  compounds,  however,  rapid  solidification  processing  is 
required.  Additional  segregants  such  as  B  can  be  deliberately  added  to  enhance  intergranular  cohesion,  and  the 
precipitation  rate  for  the  secondary  hardening  reaction  can  be  increased  to  limit  the  time  at  temperature  for  harmfu 
grain  boundary  segregation. 

Pgsign 

The  first  exercise  in  the  design  of  this  steel  is  to  estimate  the  core  and  case  c^bon  levels  required  for  the 
desired  hardness.  This  is  done  by  fitting  data  for  existing  secondary  hardening  Ni  -  Co  steels  and  extrapolating 
a  linear  fashion  to  the  desired  strength.  It  is  estimated  that  a  core  carbon  content  of  0.25  wt%  and  a  case  carbon 
content  of  0.55  wt%  is  needed  to  provide  the  desired  core  and  case  hardness  in  this  Ni  -  Co  steel.  .  ,  .  . 

The  next  step  is  to  determine  the  matrix  composition  of  the  Fe,  Ni,  and  Co.  In  order  to  produce  the  desi^ 
lath  martensite  morphology,  an  Ms  temperature  of  200  C  or  above  is  requ^ed.  Using  a  kinehc  mo^\  for  the 
compositional  dependence  of  the  Ms  temperature,  the  variation  with  Ni  and  Co  content  is  dete^ed.  This  resul 
is  illustrated  in  Figure  7  for  the  case  carbon  content  using  a  prelimmary  composition  of  the  caibide  formers  equal 


Fe-Co-Ni-5Cr-0.5Mo-0.06V-C 


Figure  7  -  Ms  temperature  calculations  for  the  case  composition  of  the  secondary-hardening 

gear  steel(Ref  [21]). 


to  5  wt%  Cr,  0.5  wt%  Mo,  and  0.06  wt%  V.  Since  the  case  has  a  higher  carbon  content  th^  the  core,  we  can  be 
assured  the  core  will  possess  a  higher  Ms  temperature  than  the  case.  In  Figure  7  the  Co  and  Ni  content  r^uired  to 
fix  the  Ms  temperature  at  200  C  is  indicated.  Since  a  high  Ni  content  is  desired  to  avoid  cleavage  fracture,  a 
composition  containing  25  wt%  Co  is  chosen.  This  allows  the  highest  possible  Ni  content,  approximate  y  3. 
wt%,  to  be  used.  These  calculations  are  later  repeated  for  consistency  when  the  composition  of  the  carbide 

formers  is  further  refined.  .  i-  j  -in.  .  *  i  •• 

To  define  the  composition  of  the  carbide  formers  a  number  of  design  constraints  are  applied.  The  total  am^nt 

of  carbide  formers  in  the  alloy  must  be  greater  than  that  r^uired  to  consume  the  carbon  prwent  ^  die  case.  This 
lower  limit  insures  that,  at  completion,  embrittling  cementite  is  completely  converted  to  M2C  carbide.  In  order  to 
reduce  grain  boundary  segregation,  the  precipitation  rate  is  maximized.  This  allows  the  shortest  possible 
tempering  time.  The  coherent  precipitation  driving  force  is  maximized  to  provide  a  small  mucal  p^cle  size  tor 
the  M2C  and  more  efficient  strengthening.  Finally,  the  solution  temperature  is  lipted  to  1000  C.  TJis  allows  Cr, 
Mo  and  V  containing  carbides  such  as  M23C6,  M7C3,  MC  and  M6C  to  be  dissolved  at  reasonable  processing 
temperatures  while  maintaining  very  fine  scale  TiC  carbides  to  act  as  the  grain  refining  dispersion.  Calculations 
for  the  precipitation  rate  constant  similar  to  those  outlined  in  the  previous  section  indicate  low  Mo  compositions 
are  favorable,  while  driving  force  calculations  have  demonstrated  the  highly  beneficid  effect  of  higher  V  co^t^ 
The  solubility  constraints  are  presented  by  the  diagram  in  Figure  8.  Here  the  equilibrium  phase  fields  at  1000  C 
are  given  as  a  function  of  Cr  and  V  content.  The  Mo  content  is  determined  by  the  stoichiometry  requirements,  the 
matrix  composition  is  taken  from  the  earlier  calculations,  and  the  carbon  content  represents  the  case  composition. 
The  point  on  the  diagram  within  the  single  phase  FCC  field  that  maximizes  the  V  content  and  nunim^s  the  Mo 
content  represents  the  best  compromise  fulfilling  all  the  design  criteria.  This  composition  is  4.8  wt%  Cr,  U.U3 
wt%  Mo  and  0  06  wt%  V.  A  recalculation  of  the  matrix  composition  using  the  final  carbide  formers  results  m  an 
alloy  composition  of  Fe  -  25  Co  -  3.8  Ni  -  4.8  Cr  -  0.03  Mo  -  0.06  V  -  0.55(case)/0.25(core)  C  (in  wt%). 
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Figure  8 


Phase  fields  on  a  pseudo-ternary  diagram  as  a  function  of  Cr, 
Mo,  and  V  content.  The  C  content  is  0.55  wt%  and  the 
remaining  composition  is  taken  from  the  matrix  composition. 


prnfotvpe  Evaluation 

A  17  lb.  vacuum  induction  heat  of  the  above  composition  was  prep^ed  by  C^enter  '^hnology  Cor^tion 
from  high  purity  materials.  The  ingot  was  forged  at  1150  C  in  a  bar  1.25  square  by  38  long  The  Ms 
temperature  of  the  alloy  was  determined  from  dilatometery  and  found  to  agree  with  rnodel  predictions.  The 
solution  treatment  response  of  the  alloy  was  determined  from  harckess  measurements  in  the  stage  I  temijr^ 
condition  The  optimum  processing  conditions  for  the  core  matenal  was  determined  to  be  a  1050  C  1  nou 
solution  treatment  followed  by  an  oil  quench  and  a  liquid  nitrogen  deep  freeze.  After  optimal  spl^on  treatment,  a 
12  hour  temper  at  482  C  results  in  the  desired  overaged  hardness  of  55  Rc  for  the  core  matenal.  The  matend  was 
then  plasma  carburized  and  processed  using  these  parameters.  The  C  potential,  temperature  and  nme  used  m  the 
carburizing  treatment  were  determined  from  the  DICTRA  software  to  provide  the  t^get  surface  carbon  content  of 
0.55  wt%  and  a  1  mm  case  depth.  Figure  9  represents  the  hardness  profile  achieved  for  the  carbimz^  sample.  A 
surface  hardness  of  67  Rc  and  a  case  depth  of  1  mm  are  obtained.  Unfortunately,  low  grmn  bound^  cohesion 
led  to  intergranular  fracture  and  limited  toughness  in  this  alloy.  Based  on  these  results,  design  work  is  continuing 
in  order  to  preserve  the  excellent  strengthening  response  and  improve  the  toughness  in  this  class  ot  alloy. 
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Figure  9  -  Carburized  hardness  profile  of  prototype  secondary  hardening  gear  steel 

(C2)  and  convention^  gear  steel(Ref  [22]). 


CONCLUSIONS  AND  SUMMARY 

The  above  example  represents  one  area  of  investigation  among  many  at  Northwestern  University  within  the 
SRG  and  Ferrium  Division.  A  list  of  materials  design  efforts,  both  previous  and  current,  follows: 

®  Stainless  bearing  steel  for  space  shuttle  main  engine  turbopumps  (patented). 

°  Aircraft  armor  steels. 

°  Oxide  superconductors. 

°  Ferritic  stainless  superalloys. 

°  Low-cost  secondary  hardening  gear  steels. 

°  Ultrahigh-strength  austenitic  TRIP  steels. 

®  Biomedcal  Co  based  corrosion  resistant  alloys. 

°  Alloys  for  high  pressure  fuel  injectors. 


Dental  concrete/Polymer  composites. 

Polymer  based  case  hardened  gears. 

Mg  based  semi-solid  casting  alloys. 
Self-Healing  biomimetic  superalloy  composites. 


The  kev  to  successfully  addressing  these  design  challenges  rests  with  the  development  of  quanutanve  models 
asdesigf^ods^Tusi^  tothedesign  of  new  matenals.  The  atove  examples 

illustrate  both  the  care  that  must  be  exercised  in  the  development  of  quantitative  useful  niodelling,  and  e 
deliberate  manner  in  which  they  need  to  be  used.  Advances  in  computational  modelling  and  computer  power 
insure  that  computational  metallurgy  will  play  an  increasing  role  in  defining  the  discipline  of  matenals  science  and 

engineering. 
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